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Abstract. We summarize the connection between random partitions and M = 2 su- 
persymmetric gauge theories in 4 dimensions and indicate how this relation extends to 
higher dimensions. 
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1. Introduction 

1.1. Random partitions. A partition of n is a monotone sequence 

A = (Ai > As > •■ • > 0) 

of nonnegative integers with sum n. The number n is denoted |A| and called the 
size of A. A geometric object associated to a partition is its diagram; it contains 
Ai squares in the first row, A2 squares in the second row and so on. An example, 
flipped and rotated by 135° can be seen in Figure ^ Partitions naturally label 




Figure 1. The diagram of A = (10,8,7,4,4,3,2,2,1,1), flipped and rotated by 135°. 
Bullets indicate the points of 6(A). The profile of A is plotted in bold. 

many basic objects in mathematics and physics, such as e.g. conjugacy classes and 
representations of the symmetric group S{n), and very often appear as modest 
summation ranges or indices. A simple but fruitful change of perspective, which 1 
wish to stress here, is to treat sums over partitions probabilistically, that is, treat 
them as expectations of some functions of a random partition. 



*The author thanks Packard Foundation for partial financial support. 
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A survey of the theory of random partitions may be found in Of the 

several natural measures on partitions, the Plancherel measure 

Mpianch(A) = , |A|=n, (1) 

stands out as the one with deepest properties and widest applications. Here dim A 
is the dimension of the corresponding representation of S{n). This is a probability 
measure on partitions of n. It can be viewed as a distinguished discretization of 
the GUE ensemble of the random matrix theory. Namely, a measure on partitions 
can be made a point process on a lattice by placing particles in positions 

e(A) = {A, - i + i} c Z + i. 

Figure n illustrates the geometric meaning of this transformation. An important 
theme of recent research was to understand how and why for a Plancherel random 
partition of n — > oo the particles 6(A) behave like the eigenvalues of a random 
Hermitian matrix. See piBl ll4[ I31j and e.g. |^,.32j^ for a survey. 

In these notes, we consider a different problem, namely, the behavior of iQ} 
in a periodic potential, that is, additionally weighted by a multiplicative periodic 
function of the particles' positions. This leads to new phenomena and new appli- 
cations. As we will see, the partition function of OTpianch in a periodic potential is 
closely related to Nekrasov partition function from supersymmetric gauge theory. 
This relationship will be reviewed in detail in Section El and its consequences may 
be summarized as follows. 

1.2. Instanton counting. In 1994, Seiberg and Witten proposed an exact 
description of the low-energy behavior of certain supersymmetric gauge theories 
|89l I4()| . In spite of the enormous body of research that this insight has gen- 
erated, only a modest progress was made towards its gauge-theoretic derivation. 
This changed in 2002, when Nekrasov proposed in a physically meaningful 
and mathematically rigorous definition of the regularized partition function Z for 
supersymmetric gauge theories in question. 

Supersymmetry makes the gauge theory partition function the partition func- 
tion of a gas of instantons. Nekrasov's idea was to use equivariant integration with 
respect to the natural symmetry group in lieu of a long-distance cut-off for the 
instanton gas. He conjectured that as the regularization parameter e ^ 

InZ - -^T 

where the free energy T expressed by the Seiberg- Witten formula in terms of 
periods of a certain differential dS on a certain algebraic curve C. 

This conjecture was proven in 2003 by Nekrasov and the author for a list 
of gauge theories with gauge group U{r), namely, pure gauge theory, theories 
with matter fields in fundamental and adjoint representations of the gauge group, 
as well as 5-dimensional theory compactified on a circle P2j. Simultaneously, 
independently, and using completely different ideas, the formal power series version 
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of Nekrasov's conjecture was proven for the pure [/(r)-theory by Nakajima and 
Yoshioka The methods of were appHed to classical gauge groups in [30] 
and to the 6-dimensional gauge theory compactified on a torus in jllj . Another 
algebraic approach, which works for pure gauge theory with any gauge group, was 
developed by Braverman and Braverman and Etingof [H]- 

In these notes, we outline the results of [221 in the simplest, yet fundamental, 
case of pure gauge theory. As should be obvious from the title, the main idea 
is to treat the gauge theory partition function Z as the partition function of an 
ensemble of random partitions. The e ^ limit turns out to be the thermodynamic 
limit in this ensemble. What emerges in this limit is a nonrandom limit shape^ an 
example of which may be seen in Figure El This is a form of the law of large 
numbers, analogous, for example, to Wigner's semicircle law for the spectrum of 
a large random matrix. The limit shape is characterized as the unique minimizcr 
tpi, of a certain convex functional S{ip), leading to 



We solve the variational problem explicitly and the limit shape turns out to be an 
algebraic curve C is disguise. Namely, the limit shape is essentially the graph of 
the function 



where dS is the Seiberg-Witten differential. Thus all ingredients of the answer 
appear very naturally in the proof. 

Random matrix theory and philosophy had many successes in mathematics 
and physics. Here we have an example when random partitions, while structurally 
resembling random matrices, offer several advantages. First, the transformation 
into a random partition problem is geometrically natural and exact. Second, the 
discretization inherent in partitions regularizes several analytic issues. For further 
examples along these lines the reader may consult 33; . 

1.3. Higher dimensions. The translation of the gauge theory problem into 
a random partition problem is explained in Section[21 In Sectional we analyze the 
latter problem, in particular, derive and solve the variational problem for the limit 
shape. Section 0] summarizes parallel results for 3-dimensional partitions, where 
similar algebraic properties of limit shapes are now proven in great generality. 

The surprising fact that free energy J- is given in terms of periods of a hidden 
algebraic curve C is an example of mirror symmetry. A general program of inter- 
preting mirror partners as limit shapes was initiated in 35 . Known results about 
the limit shapes of periodically weighted 3-dimensional partitions, together with 
the conjectural equality of Gromov-Witten and Donaldson-Thomas theories of pro- 
jective algebraic 3-folds can be interpreted as a verification of this program for 
toric Calabi-Yau 3-folds. See [3^ for an introduction to these ideas. 

Note that something completely different is expected to happen in dimensions 
> 3, where the behavior of both random interfaces and Gromov-Witten invariants 
changes qualitatively. 



J- = min S . 




dS 
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2. The gauge theory problem 

2.1. Instantons. We begin by recalling some basic facts, see |5] for an excel- 
lent mathematical treatment and [H IIOL US] for a physical one. This will serve as 
motivation for the introduction of Nekrasov's partitions function in JHJ below. 

In gauge theories, interactions are transmitted by gauge fields, that is, unitary 
connections on appropriate vector bundles. In coordinates, these are matrix-valued 
functions Ai(x) that define covariant derivatives 

d 

= —+A^{x), A* = ~A,. 

We consider the most basic case of the trivial bundle x C over the fiat Euclidean 
space-time M'', where such coordinate description is global. 

The natural (Yang-Mills) energy functional for gauge fields is L^-norm squared 
||-F||^ of the curvature 

F = ^ [Vi, Vj] dxi A dxj . 
The path integral in quantum gauge theory then takes the form 

/ PA exp(-/3||i^||2) X ... , (2) 

conncctions/C? 

where dots stand for terms involving other fields of the theory and Q is the group 
of gauge transformations g : R'' — > U{r) acting by 

In these notes, we will restrict ourselves to pure gauge theory, which is already 
quite challenging due to the complicated form of the energy. A parallel treatment 
of certain matter fields can be found in |29| . 

Our goal is to study as function of the parameter (3 (and boundary con- 
ditions at infinity, see below). A head-on probabilistic approach to this problem 
would be to make it a theory of many interacting random matrices through a dis- 
cretization of space-time. This is a fascinating topic about which I have nothing to 
say. In a different direction, when (3^0, the minima of should dominate the 
integral. In supersymmetric gauge theory, there is a way to make such approxima- 
tion exact, thereby reducing the path integral to the following finite-dimensional 
integrals. 

Local minima of are classified by a topological invariant C2 £ Z, 

called charge, and satisfy a system of first order PDEs 

F ± = , (3) 

where * is the Hodge star operator on 2-forms on R''. With the plus sign, © 
corresponds to C2 > and is called the anti-self-duality equation. Its solutions are 
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called instantons. Minima with C2 < are obtained by reversing the orientation 
of Ml 

The ASD equations lj2Jl are confornially invariant and can be transported to 
a punctured 4-sphere S"* = R** U {00} via stereographic projection. From the 
removable singularities theorem of Uhlenbeck it follows that any instanton on 
extends, after a gauge transformation, to an instanton on S^. Thus we can talk 
about the value of an instanton at infinity. 

Let Qo be the group of maps g : S'^ ^ U{r) such that g{oo) = 1. Modulo 
Qq, instantons on S*^ with C2 = n are parametrized by a smooth manifold A^(r, n) 
of real dimension 4rn. Naively, one would like the contribution from charge n 
instantons to ||2Jl to be the volume of M{r,n) in a natural symplectic structure. 
However, A^(r, n) is noncompact (and its volume is infinite) for two following 
reasons. 

Approximately, an element of Ai (r, n) can be imagined as a nonlinear superpo- 
sition of n instantons of charge 1. Some of those may become point-like, i.e. their 
curvature may concentrate in a (5-function spike, while others may wander off to 
infinity. A partial compactification of A4(r, n), constructed by Uhlenbeck, which 
replaces point-like instanton by just points of K^, takes care of the first problem 
but not the second. Nekrasov's idea was to use equivariant integration to regularize 
the instanton contributions. 



2.2. Equivariant regularization. The group 

K = SU{2) X SU{r) 

acts on A4 (r, n) by rotations of R'' = and constant gauge transformation, 
respectively. Our plan is to use this action for regularization. Let's start with the 
following simplest example: suppose we want to regularize the volume of M^. A 
gentle way to do it is to introduce a Gaussian well 

e-i''(^'+y')dxdy , m>0 (4) 

and thus an effective cut-off at the scale. Note that the Hamiltonian flow 

on R^ generated hy H — ^{x'^ + y^) with respect to the standard symplectic form 
(jj = dx A dy is rotation about the origin with angular velocity one. This makes (0J 
a simplest instance of the Atiyah-Bott-Duistermaat-Heckman equivariant localiza- 
tion formula 0. We will use localization in the following complex form. 

Let r = C* act on a complex manifold X with isolated fixed points X'^. 
Suppose that the action of C/(l) C T is generated by a Hamiltonian H with respect 
to a symplectic form uj. Then 



-2TvtH{x) 



where t should be viewed as an element of Lie(r) = C, an so it acts in the complex 
tangent space TxX to a fixed point x € X. While ^ is normally stated for compact 
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manifolds X, example (0} shows that with care it can work for noncompact ones, 
too. Scaling both lo and H to zero, we get from ((SJ a formal expression 

which does not depends on the symplectic form and vanishes if X is compact. 
A theorem of Donaldson identifies instantons with holomorphic bundles on 
— and thus gives a complex description of M {r, n) . Concretely, (r, n) is 

the moduli space of rank r holomorphic bundles f CP^ with given 2nd Chern 

class C2{£) = n and a given trivialization along the line 

ioo = \ 

at infinity. Note that existence of such trivialization implies that ci{£) = 0. A 
similar but larger moduli space Ai{r,n) of torsion-free sheaves, see e.g. [121 124| . is 
a smooth partial compactification of Ai (r, n) . 
The complexification of K 

Kc = SL{2) X SL{r) 

acts on A4{r,n) by operating on and changing the trivialization at infinity. 
Equivariant localization with respect to a general t G Lie{K) 

t ~ (diag(— ie, is), diag(iai, . . . , iar)) (7) 

combines the two following effects. First, it introduces a spatial cut-off parameter e 
as in Q). Second, it introduces dependence on the instanton's behavior at infinity 
through the parameters a^. While the first factor in K works to shepherd run- away 
instantons back to the origin, the second works to break the gauge invariance at 
infinity. In supersymmetric gauge theories, the parameters correspond to the 
vacuum expectation of the Higgs field and thus are responsible for masses of gauge 
bosons. In short, they are live physical parameters. 

2.3. Nekrasov partition function. We are now ready to introduce our 
main object of study, the partition function of the pure {J\f — 2 supersymmetric) 
U (r) gauge theory: 

Z(e;ai,...,a,;A) = Zpe,t 5]A2™ / 1, (8) 

n>0 -'Mir.n) 

where the integral is defined by © applied to (|7|), 

A = exp(-47r^/3/r) , 

and .Zport is a certain perturbative factor to be discussed below. The series in ^ 
is denoted .^inst- Because of factorials in denominators, see Zinst converges 

whenever we avoid zero denominators, that is, on the complement of 



Oi — Oj = mod s 



(9) 
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In essence, these factorials are there because the instantons are unordered. Also 
note that Z is an even function of e and a symmetric function of the a^'s. 
Since by our regularization rule 

volR*= / 1 = 4, 



we may expect that as e ^ 

lnZ(£;a;A) - -\T{a;k), 

where J- is the free energy. At first, the poles @ of Z^^st, which are getting denser 
and denser, may look like a problem. Indeed, poles of multiplicity 0{e^^) may 
affect the free energy, but it is a question of competition with the other terms in 
■^instj which the pole-free terms win if \ai — aj\ S> 0. As a result, either by passing 
to a subsequence of e, or by restricting summation in Zjnst to the relevant pole-free 
terms, we obtain a limit 

•^inst = — linie'^^inst , 

which is analytic and monotone far enough from the walls of the Weyl chambers. 
Recall that Weyl chambers for SU (r) are the r! cones obtained from 



C+ = |ai > 02 > • • • > flr, ^ tti = j 



by permuting the coordinates. As \ai — aj\ get small, poles do complicate the 
asymptotics. This is the origin of cuts in the analytic function J- (a), ai G C. 

Nekrasov conjectured in [2H| that the free energy T is the Seiberg- Witten pre- 
potential, first obtained in |39[ I4U| through entirely different considerations. It is 
defined in terms of a certain family of algebraic curves. 

2.4. Seiberg- Witten geometry. In the affine space of complex polyno- 
mials of the form P{z) — + 0{z^~'^) consider the open set U of polynomials 
such that 

P{z) = ±2K' (10) 

has 2r distinct roots. Over U, we have a g-dimensional family of complex algebraic 
curves C of genus g = r — \ defined by 

[w+^ =P{z), PeU. (11) 

The curve Hll(l is compactified by adding two points dC {w = 0, oo}. 

Let M C U be the set of P{z) for which all roots of (|10|) arc real. The corre- 
sponding curves C are called maximal and play a special role, see e.g. |4T) . They 
arise, for example, as spectral curves of a periodic Toda chain |42j . A maximal 
curve C has r real ovals, as illustrated in FigureEl Note that for z G M, w is either 
real or lies on the unit circle \'w\ = 1. 
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The intervals p-i([-2A'', 2A'']) C M on which \w\ = 1 are cahed bands. The 
intervals between the bands are called gaps. The smaller (in absolute value) root 
w of the equation can be unambiguously defined for z € C \ {bands}. On the 
corresponding sheet of the Riemann surface of w, we define cycles 

a,GH,iC^dC), l3,eHiiC,dC), i = l,...,r (12) 

as illustrated in Figure O where dotted line means that (3i continues on the other 
sheet. Note that ai n f3j — Sij and that 

a, = -ai , l3i = Pi, (13) 

where bar stands for complex conjugation. The ovals in Figure |5] represent the 
cycles ai and Pi — Pi+i- 



Pi 




: P2 
P1-P2 

— . ^ ---- 




/33 


k 










O-l 0L2 



Figure 3. Cycles f3i go from w — 00 to w = 0. Bold segments indicate bands. 

The Seiberg-Witten differential 

1 dw r , „ , , 

dS=—z — = ±— (1 + (z-^)) dz 
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is holomorphic except for a second order pole (without residue) at dC. Its deriva- 
tives with respect to P S U are, therefore, holomorphic differentials on C. In fact, 
this gives 

TpU = holomorphic diff. on C . 
Nondegeneracy of periods implies the functions 

dS, = 0, (14) 

which are real on M by (|13|l . are local coordinates on U, as are 

ar-ar_,,'^2m [ dS , J] a,^ - . (15) 

Further, there exists a function ^(a; A), which is real and convex on M, such that 

Indeed, the Hessian of J-' equals (— 27ri) times the period matrix of C, hence sym- 
metric (and positive definite on M). The function J-' is called the Seiberg-Witten 
prepotential. Note that J- is multivalued on U and, in fact, its monodromy played 
a key role in the argument of Seiberg and Witten. By contrast, M is simply- 
connected, indeed 

a"" :M^C+ 

is a diffeomorphism, see e.g. ^H] for a more general result. Note that the peri- 
ods (|15l) are the areas enclosed by the images of real ovals of C under {zjw) i-^ 
{z, In \w\). A similar geometric interpretation of the a^'s will be given in l|4U|l below. 
In particular, the range 

A = a(M) 

of the coordinates 114|) is a proper subset of C_ = — C+. At infinity of M, we have 

tti {roots of P} , Oi ^ 02 ^ • • • <C Or . 

2.5. Main result. We have now defined all necessary ingredients to confirm 
Nekrasov's conjecture in the following strong form: 

Theorem 1 ([211). For a e A, 

- \im e^\nZ{e; a; A) ^ J'ia; A) , (17) 

where T is the Seiberg- Witten prepotential (|16() . 

At the boundary of A, free energy has a singularity of the form 

where dots denote analytic terms. This singularity is one of the main physical 
features of the Seiberg-Witten theory. 

In broad strokes, the logic of the proof was explained in the Introduction. We 
now proceed with the details. 
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3. The random partition problem 

3.1. Fixed points contributions. A rank 1 torsion-free sheaf on is a 
fancy name to call an ideal / of C[a;, y]. Any partition A defines one by 

h = {x^\x^- y,x^' y^,...) (lC[x,y]. 

It is easy to see that all torus- fixed points of A4 {r, n) have the form 

r 

^? = 0^A<^)' El^^'^l^"' (18) 

fc=l 

where A^*^-* is an r-tuple of partitions. Our goal now is to compute the character 
of the torus action in the tangent space to the fixed point H18|l and thus the 
contribution of ^ to the sum in ©. 

By construction of M{r,n), its tangent space at ^ equals Extp2 (J', Loo))- 
From the vanishing of the other Ext-groups we conclude 

tre*|Ext;,(5,5(^Loo)) = -^oe-W - -^sW , (19) 
where X^{t) is the character 

-^sW =tre*|^^^(;j,5) 
of the infinite-dimensional virtual representation 

Xc- id, d) - Ext° 2 (y, - Ext ^2 (5, + Ext^ 2 . 
Any graded free resolution of 5 gives 

X^it) = \Gs{t)\\ ieLie(X), 
where G^si^) a factor, the character of ^ itself 

GA(i),...,A(^)W = (e-*^/'-e"/')tre*|s 

r oo 

= ^e-^5]exp(*£(Af -j + i)) . (20) 

k=l j=l 

It is also a natural generating function of the r-tuple A'^'^-' . 

Note that the weight of any ^ is real and positive, being a product of purely 
imaginary numbers in conjugate pairs. 

3.2. Perturbative factor. In the spirit of the original uncompactified gauge 
theory problem on R**, we would like to drop the first term in p9|) and declare its 
contribution canceled by Zpcrt- In view of (|20|l . this requires a regularization of 
the following product 

r oo 
k,k' — l j^j' — l 
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A natural regularization is provided by Barnes' double F- function (|21|l . see e.g. 
PS]. For ci, C2 e M and > 0, define 



1 f°° dt 

r(s) Jo t 



mJo t n(i-e-^'*)' 

This has a meromorphic continuation in s with poles at s = 1, 2. Define 



(21) 



Through the difference equation 

wT2iw) r2iw + Ci + C2) = r2(u; + Ci) r2(u; + C2) (22) 

it extends to a meromorphic function of w. We define 



■^port — JJ^ r 



2 

k.k' 



i{ak - ak') 
A 



A'~ 



-1 



(23) 



where r2 is analytically continued to imaginary arguments using 

1 

T2iMw\Mc,~Mc) = M2^"i2 r2 (w I c,-c) , M ^ (-cx),0]. 

The scaling by A is introduced in (|23|l to make (|SJ homogeneous of degree in a, 
e, and A. Note also 

r2(o| 1,-1) = 

Our renormalization rule 123() fits nicely with the following transformation of the 
partition function Z. 

3.3. Dual partition function. For r = 1, the weight of I\ in ||SJ) equals 



A^" det'h 



1 /A 



OTpianch(A) , (24) 



where OTlpianch is the Planchcrel measure ^ and the prefactor is the Poisson weight 
with parameter h? /e^. For r > 1, we will transform Z into the partition function 
(f^ of the Plancherel measure in a periodic potential with period r. 

Let a function ^ : Z + ^ — > M be periodic with period r and mean 0. The energy 
S(A) of the configuration (3(A) in the potential ^ is defined by Abel's rule 

S(A)= 5] e(^) = lim ^ e(^)e^^. 

Grouping the points of 6(A) modulo r uniquely determines an r-tuple of partitions 
A'-*'\ known as r-quoticnts of A, and shifts Sfe G Q such that 



6(A) = U r (6 (aC')) + ,sfe) (25) 



fc=i 
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and 



rs = p mod r Zq , P = ( ^-j^ j ■ • • j 



where Zq denotes vectors with zero sum. It foUows from (|25f) that 

Ga (e/r) = Ga(1),...,a('-) (e; es) . (26) 
Letting e 27rifc, fc = l,...,r— 1, in (|26() gives 

S(A)-(s,e)-E^'^»' (27) 
while the e hmit in H2t)|) yields 

w = -(E|^'"|+E|) + ^^ 

Using these formulas and the difference equation H22|l . we compute 

Z^(£;a,...,e.;A) ^ exp^^] Z(r£;a;A) (28) 



E 



/dim Ay /^S(A) 



exp ^ . (29) 



We call H28f) the dual partition Junction. By (|29|l . it equals the partition function 
of a periodically weighted Plancherel measure on partitions. 

While it will play no role in what follows, it may be mentioned here that 
is a very interesting object to study not asymptotically but exactly. For example, 
Toda equation for In Z ^ may be found in Section 5 of ■ 

3.4. Dual free energy. Define the dual free energy by 

J^^(C;A) = - lime^lnZ^ . (30) 

Since (|28|l is a Riemann sum for Laplace transform, we may expect that 

^^(^;A) = minl.F(a;A)-i(e,a) (31) 

that is, up to normalization, \s the Legendre transform of !F. This is because 
the asymptotics of Laplace transform is determined by one point — the maximum. 
Our plan is apply to same logic to the infinite-dimensional sum H29II . namely, to 
show that its e ^ asymptotics is determined by a single term, the limit shape. 

The law of large numbers, a basic principle of probability, implies that on a large 
scale most random system are deterministic: solids have definite shape, fluids obey 
the laws of hydrodynamics, etc. Only magnification reveals the full randomness of 
nature. 
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In the case at hand, the weight of a partition A in (|29|) . normahzed by the whole 
sum, defines a probability measure on the set of partitions. This measure depends 
on a parameter e and as £ —> it clearly favors partitions of larger and larger size. 
In fact, the expected size of A grows as e~^. We thus expect the diagram of A, 
scaled by e in both directions, to satisfy a law of large numbers, namely, to have a 
nonrandom limit shape. By definition, this limit shape will dominate the leading 
e ^ asymptotics of . In absence of the periodic potential S, such analysis is 
a classical result of Logan-Shepp and Vershik-Kerov |21l 1431 Ej . 

Note that the maximum in 1)31(1 is over all of a, including the problematic 
region where \ai — aj\ get small. However, this region does not contribute to 
as the convexity of free energy is lost there. We will see this reflected in the 
following properties of JF^; it is strictly concave, analytic in the interior of the 
Weyl chambers, and singular along the chambers' walls. 

3.5. Variational problem for the limit shape. The profile of a par- 
tition A is, by definition, the piecewise linear function plotted in bold in Figure ^ 
Let be the profile of A scaled by e in both directions. The map ^\ embeds 
partitions into the convex set ^ of functions on M with Lipschitz constant 1 and 

V'l = j \^{x) — \x\\dx < cxo . 

The Lipschitz condition implies 

1101 -02||C< 1101- ^2|li(', 01,0^2 e*, 

and so \& is complete and separable in the L^-metric. Some function of a partition 
have a natural continuous extension to for example 

Ga(^) - ,../2 ^,-.e/2 (l - ^/ (^a(-) - 

while others, specifically the ones appearing in H29|l . do not. An adequate language 
for dealing with this is the following. 

Let /(A) > be a function on partitions depending on the parameter e. We 
say that it satisfies a large deviation principle with action (rate) functional Sf{4>) 
if for any set A C '4' 



-lime^ In ^ /(A) C 



inf Sf, inf 5f 



cMU{+cx)}, (32) 



where lim denotes all limit points, A° and A stand for the interior and closure of 
A, respectively. 

For the Plancherel weight (|24|l . Logan-Shepp and Vershik-Kerov proved a large 
deviation principle with action 

'5pi(0)-^/ (l-h0'(a;))(l-0'(y)) In^-^dxdy. (33) 

^ J x<y 
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Note that in this case the sum in (|32|l may be replaced by maximum because the 
number of partitions of n grows subexponentiaUy in n. In other words, there is no 
entropic contribution in (|33ll . 

The periodic potential S(A) produces a surface tension addition to the total 
action S 

where ct is a convex piecewise-linear function of the kind plotted in Figure 01 It is 
linear on segments of length 2/r with slopes {^i}, in increasing order. The form of 
5surf is easy to deduce directly; it can also be seen as e.g. the most degenerate case 
of the surface tension formula from Section 14.21 The singularities in the surface 
tension a are responsible for facets, that is, linear pieces, in the minimizer, see Fig- 
ure|H| The slopes of these facets are precisely the points where a' is discontinuous^. 
Note that cr and hence 5 is a symmetric function of the ^.^'s. 




Figure 4. The surface tension a for r — A and ? = { — |,— fif,!} 

The functional S is strictly convex and its sublevel sets <S~^((— oo, c]) are com- 
pact, which can be seen by rewriting iSpi in terms of the Sobolev H^^'^ norm, see 
[211 144|. Therefore it has a unique minimum ■0* — the limit shape. The large 
deviation principle and the definition of the dual free energy and H^Ofl together 
imply 

T^[^-K)=S{^.). (34) 
Our business, therefore, is to find this minimizer if^i,. 

3.6. The minimizer. By convexity, a local minimum of S is automatically 
a global one. Since a has one-sided derivatives, a local minimum can be charac- 
terized by nonnegativity of all directional derivatives. This leads to the following 
complementary slackness conditions for the convolution of ■0* (2^) with the kernel 

\x\ /"^ 

L(a;)=a;ln— x— / In 

A 7o 

^ There are many advantages in viewing random partitions as 2-dimensional slices of random 
3-dimensional objects discussed in Section |3 Prom the probabihty viewpoint, this links random 
partitions with rather realistic models of crystalline surfaces with local interaction, enriching 
both techniques and intuition. In particular, coexistence of facets and curved regions in our 
limit shapes is the same phenomenon as observed in natural crystals. From the gauge theory 
viewpoint, it is also very natural, especially in the context of 5-dimensional theory on X , 
which corresponds to the X-theory of the instanton moduli spaces. 



dy. 
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There exists a constant cq, which is the Lagrange multipher from the constraint 
J Sip' = 0, such that 



-1 



(35) 



L*V"(a;)+co = , 

L * ij':{x) + Co e , iPUx) = -l + f 

where to simplify notation we assumed that 

C e C_ , Co = -oo , Cr+l +00 ■ 

Recall that C_ denotes the negative Weyl chamber. 

The function will turn out to be nonnegative and supported on a union of r 
intervals, which are precisely the bands of Section [T^ The gaps will produce the 
facets in the limit shape. 

It is elementary to see that for a maximal curve ()ll|l the map 



2 2 A 

$(z) = 1 + lnw = H ln- + 0(z~^) 

TTir ni z 



(36) 



where w is the smaller root of Hll|l , defines a conformal map of the upper half-plane 
to a slit half-strip 

A c > 0,|5R^| < 1} 

as in Figure El The slits in A go along 

Kz = -l + 2i/r, i = l,...,r-l, 

and their lengths are, essentially, the critical values of the polynomial P{z). The 
bands and gaps are preimages of the horizontal and vertical segments of (3A, re- 
spectively. 



Figure 5. Conformal map defined by a maximal curve 
We claim that 



(37) 



where the polynomial P{z) is determined by the relation (|39(l below. The equations 
H35|) are verified for H37I) as follows. Since ^'{z) = 0{z~^), z oo, we have the 
Hilbert transform relation 



P.V. - * 

X 
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Integrating it once and using to fix the integration constant, we get 

(L*5R$')' = 71-3$. 

Therefore, the function L * 3?$' is constant on the bands and strictly increasing on 
the gaps, hence (pTzl satisfies with 

(38) 



- = / dx 

J ith gap 

Integrating by parts and using definitions from Section lOI gives 



e = — , 



(39) 



thus every limit shape ip^, comes from a maximal curve. For example, the limit 
shape corresponding to the curve from Figure [3 is plotted in Figure |^ Note also 




Figure 6. Limit shape corresponding to the curve from Fig.|21 Thin segments are facets, 
that for C e M, we have 



2 (^.-1 - h) , 



(40) 



where li is the intercept of the ith facet of the limit shape. In particular, A C C_. 

For given ^ G C_, consider the distribution of the r-quotients A*-*-* of the parti- 
tion A, as defined in Section For the shifts Sk in (|25|) we have using (|27|l 



es 



in probability. Observe that 



(A) 



since the other term, containing ^'Z'*, vanishes by the definition of a maximum. 
Definitions and integration by parts yield 

d d 



By (|26f) . this means that the resulting sum over the r-quotients A^*'' is the original 
partition function Z with parameters a G A. This concludes the proof. 
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4. The next dimension 
4.1. Stepped surfaces. 



An obvious 3-dimensional generalization of a partition, 
also known as a plane partition can be seen on the right. 
More generally, we consider stepped surfaces, that is, 
continuous surfaces glued out of sides of a unit cube, 
spanning a given polygonal contour in M"^, and project- 
ing 1-to-l in the (1,1,1) direction, see Figure |3 Note 
that stepped surfaces minimize the surface area for given 
boundary conditions, hence can be viewed as zero tem- 
perature limit of the interface in the 3D Ising model. 




The most natural measure on stepped surfaces is the uniform one with given 
boundary conditions, possibly conditioned on the volume enclosed. It induces 
Plancherel-like measures on 2-dimensional slices. Stepped surfaces are in a natural 
bijection with fully packed dimers on the hexagonal lattice and Kasteleyn theory 
of planar dimers |16| forms the basis of most subsequent developments. 

The following law of large numbers for stepped surfaces was proven in . Let 
Cn be a sequence of boundary contours such that each C„ can be spanned by at 
least one stepped surface. Suppose that n~^C„ converge to a given curve C C M'^. 
Then, scaled by n~^, uniform measures on stepped surfaces spanning C„ converge 
to the (5-measure on a single Lipschitz surface spanning C — the limit shape. This 
limit shape formation is clearly visible in Figure [T] 




Figure 7. A limit shape simulation. The frozen boundary is the inscribed cardioid. 
The limit shape is the unique minimizer of the following functional. Let the 
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surface be parameterized by X3 = h{x3 — xi,X3 — X2), where h is a Lipschitz 
function with gradient in the triangle A with vertices (0, 0), (0, 1), (1, 0). Let be 
the planar region enclosed by the projection of C in the (1,1,1) direction. We will 
use {x, y) — {xz — xi, X2 — xi) as coordinates on Q. The limit shape is the unique 
minimizcr of 

Sstcp{h) = / astcpi'^h)dxdy , (41) 

where, in the language of '20' , the surface tension (Tstcp is the Legendre dual of the 
Ronkin function of the straight line 

z + w = 1. (42) 

We recall that for a plane curve P(z, w) = 0, its Ronkin function is defined by 

„^ X 1 f f , dz dw , , „, 

{2Tny JJ ^^-"^ ' z w 

\w\—e'^ 

The gradient VR always takes values in the Newton polygon A(P) of the polyno- 
mial P, so A(P) is naturally the domain of the Legendre transform i?^. For the 
straight line as above, the Newton polygon is evidently the triangle A. 

The surface tension (Tstep is singular and not strictly convex at the boundary of 
A, which leads to formation of facets and edges in the limit shape (which can be 
clearly seen in Figure [7J. This models facet formation in natural interfaces, e.g. 
crystalline surfaces, and is the most interesting aspect of the model. Note that 
facets are completely ordered (or frozen). The boundary between the ordered and 
disordered (or liquid) regions is known as the frozen boundary. 

The following transformation of the Euler-Lagrange equation for 141|) found 
in 2^ greatly facilitates the study of the facet formation. Namely, in the liquid 
region we have 

V/i — — (arg w, — arg z) , (44) 

TT 

where the functions z and w solve the differential equation 

^ + ^ = c (45) 

z w 

and the algebraic equation H42|) . Here c is the Lagrange multiplier for the volume 
constraint J^h = const, the unconstrained case is c = 0. At the boundary of 
the liquid region, z and w become real and the \7h starts to point in one of the 
coordinate directions. 

The first-order quasilinear equation 145() is, essentially, the complex Burgers 
equation z^ ~ zzy and, in particular, it can be solved by complex characteristics 
as follows. There exists an analytic function Q{z, w) such that 

g(e-'=^z,e-"2'u;) = 0. (46) 

In other words, z(x, y) can be found by solving (|42(l and (|46(l . In spirit, this is very 
close to Weierstrafi parametrization of minimal surfaces in terms of analytic data. 
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Frozen boundary can only develop if Q is real, in which case the roots (z, w) 
and (z, w) of H46I) coincide at the frozen boundary. At a smooth point of the frozen 
boundary, the multiplicity of this root will be exactly two, hence Vft. has a square- 
root singularity there. As a result, the limit shape has an x^/^ singularity at the 
generic point of the frozen boundary, thus recovering the well-known Pokrovsky- 
Talapov law |37| in this situation. At special points of the frozen boundary, triple 
solutions of (|46|l occur, leading to a cusp singularity. One such point can be seen 
in Figured 

Remarkably, for a dense set of boundary condition the function Q is, in fact, 
a polynomial. Consequently, the frozen boundary takes the form i?(e^^,e'^^) = 0, 
where R is the polynomial defining the planar dual of the curve Q — Q. This 
allows to use powerful tools of algebraic geometry to study the singularities of the 
solutions, see ^^l- The precise result proven there is 

Theorem 2 f|19|'l. Suppose the boundary contour C is a connected polygon with 
3fc sides in coordinate directions (cyclically repeated) which can be spanned by a 
Lipschitz function with gradient in A. Then Q = is an algebraic curve of degree 
k and genus zero. 

For example, for the boundary contour in Figure H we have fc = 3 (one of the 
boundary edges there has zero length) and hence R is the dual of a degree 3 genus 
curve — a cardioid. The procedure of determining Q from the boundary conditions 
is effective and can be turned into a practical numeric homotopy procedure, see 
|19| . Higher genus frozen boundaries occur for multiply-connected domains, in 
fact, the genus of Q equals the genus of the liquid region. 

Of course, for a probabilist, the law of large numbers is only the beginning and 
the questions about CLT corrections to the limit shape and local statistics of the 
surface in various regions of the limit shape follow immediately. Conjecturally, 
the limit shape controls the answers to all these questions. For example, the 
function e~'^^z defines a complex structure on the liquid region and, conjecturally, 
the Gaussian correction to the limit shape is given by the massless free field in 
the corresponding conformal structure. In the absence of frozen boundaries and 
without the volume constraint, this is proven in ^7]. See e.g. [T51 [T71 OUllS^l ISl] 
for an introduction to the local statistics questions. 

4.2. Periodic weights. Having discussed periodically weighted Plancherel 
measure and a 3-dimensional analog of the Plancherel measure, we now turn to 
periodically weighted stepped surfaces. This is very natural if stepped surfaces are 
interpreted as crystalline interfaces. Periodic weights are introduced as follows: 
we weight each square by a periodic function of 2:3 — xi and X2 — xi (with some 
integer period M). 

The role previously played by the straight line (|42|l is now played by a certain 
higher degree curve P{z,w) = 0, the spectral curve of the corresponding periodic 
Kasteleyn operator. In particular, the surface tension dstep is now replaced by the 
Legendre dual of the Ronkin function of P, see |201- We have 



deg P = M 
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and the coefRcients of P depend polynomially on the weights. 

The main result of |2()j . known as maximality, says that for real and positive 
weights the curve P is always a real algebraic curve of a very special kind, namely, 
a Harnack curve, see |23| . Conversely, as shown in |18j , all Harnack curves arise 
in this way. 



Figure 8. The (curved part of the) Ronkin function of a genus 1 curve. Its projection to 
the plane is the amoeba. 

Harnack curves are, in some sense, the best possible real curves; their many 
remarkable properties are discussed in One of several equivalent definitions 

of a Harnack curve is that the map 

i-^ (log|z|,log|w|) (47) 

from P{z,w) to is 1-to-l on the real locus of P and 2-to-l over the rest. 
The image of P = under (|47|l is known as the amoeba of P. Note from H43|l 
that the gradient Vi? of the Ronkin function of P is nonconstant precisely for 
(x, y) in the amoeba of P. In other words, the Ronkin function has a facet (that 
is, a linear piece) over every component of the amoeba complement. The 2-to- 
1 property implies that the number of compact facets of Ronkin function equals 
the (geometric) genus of the curve P. Each of these facets translates into the 
singularity of the surface tension and, hence, into facets with the same slope in 
limit shapes. 

By Wulff's theorem, the Ronkin function itself is a minimizer, corresponding 
to its own ("crystal corner") boundary conditions. An example of the Ronkin 
function of a genus 1 Harnack curve can be seen in Figure |H1 

Maximality implies persistence of facets, namely, for fixed period M , there 
will be (*^2~^) compact facets of the Ronkin function and (*^2"^) corresponding 
singularities of the surface tension, except on a codimension 2 subvariety of the 
space of weights. It also implies e.g. the following universality of height fluctuations 
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in the liquid region 

Var(/i(a) - h{b)) - - In \\a - b\\ , \\a - 6|1 ^ oo . 

TT 

Remarkably, formulas (|44l) . H45I) . and (|46|l need no modifications for periodic 
weights. Replacing by P{z, w) — is the only change required, see [T^ . 

From our experience with periodically weighted Plancherel measure, it is natu- 
ral to expect that, for some special boundary conditions, the partition function of 
periodically weighted stepped surfaces will encode valuable physical information. 
A natural choice of "special boundary conditions" are the those of a crystal corner, 
when we require the surface to be asymptotic to given planes at infinity, as in Fig- 
ure |S1 For convergence of the partition function, one introduces a fugacity factor 
(7™\ where the missing volume is measured with respect to the "full corner". 




Figure 9. Two 3D partitions connected at an angle through an infinite leg. 

I hope that further study will reveal many special properties of such crystal 
corner partition functions. Their extremely degenerate limits have been identified 
with all-genera, all-degree generating functions for Donaldson-Thomas invariants 
of toric Calabi-Yau threefolds. Namely, as the periodic weights become extreme, 
all limit shapes, and the Ronkin function in particular, degenerate to piecewise 
linear functions. This is known as the tropical limit. The only remaining features 
of limit shapes are the edges and the triple points, where 2 and 3 facets meet, 
respectively. In this tropical limit, the partition function becomes the partition 
ftinction of ordinary, unweighted, 3D partitions located at triple points. These 3D 
partitions may have infinite legs along the edges, as in Figure El and through these 
legs they interact with their neighbors. This description precisely matches the 
localization formula for Donaldson-Thomas invariants of the toric threefold whose 
toric polyhedron is given by the piecewise linear limit shape, see |22| . 

Donaldson- Thomas theory of any 3-fold has been conjectured to be equivalent, 
in a nontrivial way, to the Gromov-Witten theory of the same 3-fold in 22 . For 
the toric Calabi-Yau 3-folds, this specializes to the earlier topological vertex con- 
jecture of 0. It is impossible to adequately review this subject here, see (SEl for 
an introduction. This is also related to the supersymmetric gauge theories con- 
sidered in Section[21 or rather their 5-dimensional generalizations, via a procedure 
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called geometric engineering of gauge theories. See for example and references 
therein. 

I find such close and unexpected interaction between rather basic statistical 
models and instantons in supersymmetric gauge and string theories very exciting 
and promising. The field is still full of wide open questions and, in my opinion, it 
is also full of new phenomena waiting to be discovered. 
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